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ABSTRACT 

We compare new maps of the hot gas, dark matter, and galaxies for IE 0657-56, a cluster with a rare, high- 
velocity merger occurring nearly in the plane of the sky. The X-ray observations reveal a bullet-like gas sub- 
cluster just exiting the collision site. A prominent bow shock gives an estimate of the subcluster velocity, 4500 
kms"', which lies mostly in the plane of the sky. The optical image shows that the gas lags behind the subclus- 
ter galaxies. The weak-lensing mass map reveals a dark matter clump lying ahead of the collisional gas bullet, 
but coincident with the effectively collisionless galaxies. From these observations, one can directly estimate the 
cross-section of the dark matter self-interaction. That the dark matter is not fluid-like is seen directly in the X- 
ray - lensing mass overlay; more quantitative limits can be derived from three simple independent arguments. 
The most sensitive constraint, a/m < 1 cm^g"', comes from the consistency of the subcluster mass-to-light 
ratio with the main cluster (and universal) value, which rules out a significant mass loss due to dark matter 
particle collisions. This limit excludes most of the 0.5-5 cm^g"' interval proposed to explain the flat mass 
profiles in galaxies. Our result is only an order-of-magnitude estimate which involves a number of simplifying, 
but always conservative, assumptions; stronger constraints may be derived using hydrodynamic simulations of 
this cluster. 

Subject headings: dark matter — galaxies: clusters: individual (1E0657-56) — galaxies: formation — large 
scale structure of universe 



1. INTRODUCTION 

IE 0657-56, one of the hottest and most X-ray luminous 
galaxy clusters known, was discovered by Tucker et al. 
(1995). It was first observed by Chandra in October 2000 
for 24 ks. That observation revealed a bullet-like, relatively 
cool subcluster just exiting the core of the main cluster, with 
a prominent bow shock (Markevitch et al. 2002, hereafter 
M02). A comparison of the X-ray and optical images revealed 
a galaxy subcluster just ahead of the gas "bullet", which led 
M02 to suggest that this unique system could be used to de- 
termine whether dark matter is collisional or colUsionless, if 
only one could map the mass distribution in the subcluster. 
Apart from the obvious interest for the still unknown nature 
of dark matter, the possibility of it having a nonzero self- 
interaction cross-section has far-reaching astrophysical impli- 
cations (Spergel & Steinhardt 2000; for more discussion see 
^below). 

Just such a map of the dark matter distribution in 
IE 0657-56 has recently been obtained by Clowe, Gonza- 
lez, & Markevitch (2004, hereafter C04) from weak lensing 
data. It reveals a dark matter clump coincident with the cen- 
troid of the galaxies (Fig. ^). C04 also derived M /L ra- 
tios of the main cluster and the subcluster and found them in 
agreement with each other and with other clusters' values. In 
addition, Chandra re-observed IE 0657-56 for 70 ks in July 
2002, from which a more accurate estimate of the shock Mach 
number was derived using the gas density jump at the shock, 
M = 3.2;|;Qg (all uncertainties 68%), which corresponds to a 
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shock (and bullet subcluster) velocity of Vs = 4500;!;g()o kms" 
(Markevitch et al., in prep., hereafter M04). The new X-ray 
data also further clarified the geometry of the merger. 

In this paper, we combine these new optical and X-ray data 
to constrain the self-interaction cross-section of dark matter 
particles. We use = 0.3, 0.^ = 0.7, //o=70 km s"' Mpc"', 
for which 1" = 4.42 kpc at the cluster redshift z = 0.296. 

2. COLLISIONAL CROSS-SECTION ESTIMATES 

The dark matter collisional cross-section, a, can be con- 
strained from the IE 0657-56 data by at least three indepen- 
dent methods, using simple calculations described in the sec- 
tions below. They are based on the observed gas-dark mat- 
ter offset, the high subcluster velocity, and the subcluster sur- 
vival. First, we give the main assumptions that will go into 
these calculations. 

There are two estimates of the total masses of the subclus- 
ter and the main cluster — from the galaxy velocity dispersion 
(Barrena et al. 2002) and from weak lensing (C04). Given the 
disturbed state of this system, virial or hydrostatic mass esti- 
mates (either from galaxy velocities or the gas temperature) 
can be incorrect, and we chose to use the direct weak lensing 
measurements from C04 even though their formal statistical 
accuracy is poorer The main cluster's lensing signal can be 
fit by a King mass profile p = po{l + r^/r^)"-'/^ with best-fit 
parameters po ~ 2.6 x 10"^^ gem"-' and r^ ~ 210 kpc (C04). 
These two parameters are degenerate so their individual er- 
ror bars are not meaningful; the quantity of interest to us is 
the central mass column density (approximately proportional 
to Pore), which is measured with a 16% accuracy. This mass 
profile is very close to the Barrena et al. NEW profile at all 
radii outside the core. A King profile is marginally preferred 
over an NEW profile (also acceptable statistically). 

The projected mass excess created by the subcluster is de- 
tected in the lensing data with a 3.0(7 significance. The sub- 
cluster mass signal is detected to r ~ 150-200 kpc from the 
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Fig. 1 . — (a) Overlay of the weak lensing mass contours on the optical image of 1E0657-56. Dashed contours are negative (relative to an arbitrary zero level). 
The subcluster's DM peak is coincident within uncertainties with the centroid of the galaxy concentration, (fo) Overlay of the mass contours on the X-ray image 
(only the upper 4 of those in panel a are shown for clarity). The gas bullet lags behind the DM subcluster. 



mass (or galaxy) peak, beyond which the subcluster may be 
tidally stripped (C04). The lensing-derived subcluster mass 
(M~7x IO'^Mq) is significantly higher than the Barrena et 
al. estimate, but the latter result was based on only 7 galax- 
ies and the equilibrium assumption, thus could easily be bi- 
ased. The current lensing data accuracy is not sufficient to 
derive an exact mass distribution for the subcluster, but our es- 
timates below are not particularly sensitive to it, being mostly 
determined by its overall projected mass. For the sake of 
modeling, we will adopt a King profile with rc = 70 kpc and 
po = 1 -3 X 10"^"^ g cm""', truncated at rn = 150 kpc, which ade- 
quately describes the lensing data. 

The subcluster is assumed to have passed (once) close to the 
center of the main cluster This is supported by the X-ray im- 
age (Fig.fn?), the gas temperature map (M04) and their com- 
parison with the radio halo map (Govoni et al. 2004; Liang et 
al. 2000). A cooler North-South bar in the X-ray image be- 
tween the two dark matter clumps appears to be an edge-on 
pancake-like remnant of the merged main cluster's gas core 
and the subcluster's outer atmosphere (stripped from what is 
now the gas bullet), suggesting that the subcluster has passed 
straight through the densest cluster region (M04). The line- 
of-sight velocity of the subcluster relative to the main cluster 
is about 600 kms"' (Barrena et al. 2002); combined with the 
X-ray-measured Mach number, it gives an angle of only ^8° 
between the direction of motion and the plane of the sky. The 
sharpness of the shock front also confirms that the subcluster 
is presently moving very nearly in the plane of the sky (M02). 
From all the above, it is reasonable to assume that the sub- 
cluster has passed through the core of the main cluster 

The accuracy of our qualitative cross-section estimates will 
be determined by the validity of this and other assumptions 
(given below where relevant) to a greater degree than by the 
measurement uncertainties, so below we will omit the mea- 



surement error propagation for clarity. 

2.1. The gas — dark matter offset 

The most remarkable feature in Fig.^ is a 23" offset be- 
tween the subcluster's DM centroid and the gas bullet, which 
is at least 2(T-significant (C04). C04 use this fact as a direct 
proof of dark matter existence, as opposed to modified grav- 
ity hypotheses (Milgrom 1983 and later works) in which one 
would expect the lensing mass peak to be associated with the 
gas — the dominant visible mass component. For our pur- 
poses, this offset means that the scattering depth of the dark 
matter subcluster w.r.t. collisions with the flow of dark mat- 
ter particles cannot be much greater than 1. Otherwise the 
DM subcluster would behave as a clump of fluid, experienc- 
ing stripping and drag deceleration, similar to that of the gas 
bullet (assuming the same gas mass fraction in the main clus- 
ter and the subcluster), and there would be no offset between 
the gas and dark matter The subcluster's scattering depth is 

Ts=-Ss, (1) 

m 

where a is the DM collision cross-section, m is its particle 
mass, and is the DM mass surface density of the subclus- 
ter. The surface density averaged over the face of the sub- 
cluster within r = rtr is Es ~ 0.2 gcm~^. Assuming spherical 
symmetry and requiring that Ts < 1 , we obtain 

-<5cm2g-'. (2) 
m 

The surface density toward the subcluster center is several 
times higher, so by using an average we obtain a conserva- 
tive upper limit. 

Another remarkable feature in Fig.QJ? is the coincidence of 
the subcluster's DM and galaxy centroids within their uncer- 
tainties (C04). To avoid an easily made mistake, we should 
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note here that the subcluster galaxies, although effectively 
collisionless, do not give us the position at which the sub- 
cluster would be in the absence of the DM collisions. The 
galaxies are gravitationally bound to the DM clump and will 
be dragged back by it, should it experience any deceleration 
from the DM collisions. This can be taken into account, and 
a limit on the offset between the DM and galaxy centroids 
can indeed be used to derive an independent a jm constraint; 
however, it requires more accurate centroid positions than cur- 
rently available. 

2.2. The high velocity of the subcluster 

The observed velocity of the subcluster, Vs = 4500 kms"', 
is in good agreement with the expected free-fall velocity onto 
the main cluster For our main cluster's mass profile, a small 
subcluster falling from a large distance should acquire 4400 
kms"' at core passage, decelerating to 3500 km s~' at the cur- 
rent 0.66 Mpc off-center distance of the subcluster Since 
the cluster peculiar velocities are small, such an agreement 
strongly suggests that the subcluster could not have lost much 
of its momentum to drag forces. Drag is created by puta- 
tive DM particle collisions, as well as the gravitational pull of 
the gas being stripped from the subcluster and of the tidally- 
stripped outer subcluster mass, and by dynamical friction as 
the moving subcluster disturbs the main cluster's matter dis- 
tribution. We will conservatively disregard the latter three (all 
of which are relatively small effects) and assume, for a qual- 
itative estimate, that the loss of velocity due to the DM colli- 
sions, compared to free fall, is less than 1000 kms"', accumu- 
lated along the way through the main cluster 

We will now derive the drag force on the subcluster from 
the DM particle collisions. The subcluster will be decelerated 
by collisions when its particles acquire a momentum compo- 
nent opposite to the subcluster's velocity and then transfer all 
or part of it to the whole subcluster via gravitational interac- 
tions. We assume for simplicity that (a) dark matter particles 
have no peculiar velocities, and (b) the subcluster's gravity is 
felt by the particles, and vice versa, as long as the particle is 
within a fiducial radius r' = Ir^i- of the subcluster (the result 
will not depend qualitatively on / in the range rtr-oo). We 
further assume that a particle transfers all of its momentum 
to the subcluster if its post-collision velocity is insufficient to 
escape beyond r'. Faster-scattered particles retard the cluster 
only until they reach r' . The critical velocity that a subcluster 
particle needs to reach r', mass-averaged (in the nns sense) 
over the subcluster mass profile, is V ~ 1900 kms"'. Use of 
this average value lets us make a further simplifying assump- 
tion that all collisions occur at the subcluster center We also 
conservatively disregard multiple scattering in the subcluster; 
the adequacy of this assumption will be addressed below. 




Fig. 2. — Collision of two equal-mass particles in the subcluster reference 
frame. 



An elastic collision of two equal-mass particles proceeds as 
shown in Fig.|2l In the subcluster's reference frame, particle 
2 is at rest and particle 1, in the incoming flow, collides with 
it with a velocity 



= A/v2 + y2w4800kms" 



(3) 



where Vs is the subcluster's velocity and the relatively small 
increase is because of the subcluster's gravitational pull 
within r = r' . Particle 1 scatters at an angle ai > with 
a velocity vi, while particle 2 acquires a velocity V2 at an 
angle 02 > 0. From energy and momentum conservation, 
ai + a2 = 7r/2 and the velocities are 

vi=vocosai, V2 = V()Sinai. (4) 

As a result of each collision, the subcluster acquires a net 
momentum along the vq direction, p = po + p\+ pi, where 



po = m(v.s-vo) 



(5) 



comes from the infalling particle (which carries back exactly 
-po if it passes through the subcluster without scattering), and 
Pi (i =1,2) comes from the respective outbound particle. If 
Vi < V flien 

Pi = mvi cos ai (6) 

(the particle transfers all its momentum to the subcluster), 
while for v, > V, 



Pi = m v,-(v^-y^)'''^ coso;,- 



(7) 



For our particular mass and velocity, at least one particle al- 
ways escapes beyond r' . If both escape, which occurs in col- 
lisions with Qfe < cki < f - tte (o^e — 23° for our parameters), 
the momentum loss by the subcluster as a function of the scat- 
tering angle ai is (combining eqs.|5j0|3j andO: 



/^both = mv, 



y2 



1-cosai I cos a 1 j"^'*^ 



-smai 



.2 2 

sm a] jcos a\ 



1/2 



1/2- 



(8) 



If only one particle escapes, which occurs when < ai < ae 
or f-ae < a I < ?, the net momentum loss is (combining 
.ElilElEl andgj: 



eqs. 



l-cosai cos a 



y . 2 

n — ^ sm a\ 
vi 



1/2- 



(9) 



Now we should average the above momentum loss over all 
scattering angles. For this, it is convenient to use the refer- 
ence frame of the center of mass of the colliding particles. 
In this frame, the scattering is isotropic, as long as the parti- 
cles are "slow" in the sense that mur/h ^ 1, where r is the 
linear scale of the interaction and u = v'o/2 is the collision ve- 
locity (Landau & Lifshitz 1958, §132). For example, for rigid 
spheres with radius r (for which the scattering cross-section is 
(7 = Landau & Lifshitz 1958), the scattering is isotropic 
if 



/ 2 \3/2 

/ mc \ / u \ 
1 GeV J V 2400 kms"' J 



a /m 



50cm2g 1 



1/2 



< 1 (10) 



Note that if the DM particles are WIMPs with the currently 
favored mass {mc^ > 40 GeV, Hagiwara et al. 2002), this con- 
dition is not satisfied for the values ofcr/m'^(l-10) cm^g"' 
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that we are able to constrain. Scattering with such combina- 
tions of a and m would be strongly beamed, reducing all the 
observable effects for a given a. However, there is still ample 
parameter space for isotropic scattering (e.g., if particles are 
not WIMPs), and we will confine our estimates to this simple 
case. Our conclusions will be directly comparable to other as- 
trophysical cross-section measurements which also implicitly 
assume isotropic scattering. 

In this reference frame, particle 1 scatters by an angle 9 = 
2a 1 (0 < 6* < tt) and particle 2 scatters by tt-O. Combining 
eqs. (|8l|9}, weighting with solid angle and taking into account 
the symmetries, the average momentum lost by the subcluster 
in each particle collision is (for our values of V and Vj): 



-4/'.r{/-^;;(i. 

•/ sin Oe k s 



1/2 



dx 



: 0.1 mvj. 



(11) 

When the subcluster with mass Ms (all in the form of dark 
matter for clarity) travels with a velocity v through the main 
cluster with density p^, it experiences 



M, (J 

n= pn 

m m 



(12) 



collisions per second. As a result, it loses velocity, relative to 
the free-fall velocity Vff : 



d{v-vs) 
dt 



pn 



P o- 



(13) 



where p is from eq. il 1\ . We conservatively disregard the 
higher subcluster velocity during core passage and assume 
that the loss of mass and velocity is relatively small, so p 
is roughly constant. Integration of eq. jl3t along the sub- 
cluster trajectory gives the total velocity loss. Noting that 
J Pmvdt = J Pmdl is the mass column density of the main 
cluster along the subcluster's trajectory, Sm, we obtain 



P sr. 

- Vff = 2j„ 

m m 



(14) 



Almost all of accumulates within the main cluster's 
core, at distances smaller than the current subcluster posi- 
tion. If, as we assume, the subcluster passed through the main 
cluster center, Sm ~ 0.3 gcm~^; if it missed the center by 
200-300 kpc. En, is lower by about a factor of 2. On the 
other hand, we observe the main cluster after the collision; if 
it used to have an NFW-type peak disrupted by a direct hit of 
the subcluster, then that we should use would be higher 
than the observed one. 

Requiring that v- Vff < 1000 kms"\ from eq. (I14> we get 



— < 7 cm" j 
m 



(15) 



2.3. The survival of the dark matter subcluster 



For the observed subcluster mass and velocity, the most 
likely result of a particle collision is the loss of a particle 
by the subcluster. The opposite effect, i.e., accretion of the 
main cluster's particles, is negligible because of the low mass 
and high velocity of the subcluster. We can put an upper 
limit on the integrated mass loss and thereby on the colli- 
sion cross-section. C04 have derived mass-to-light ratios for 
the subcluster within r^ = 150 kpc, M/Lg ~ 280 ±90 and 
M/Lj ~ 170 ±50. These ratios are in good agreement with 
the universal cluster values from the lensing analyses (e.g., 
Melher 1999; Dahle 2000), and a factor of 1.1 ±0.3 from the 



main cluster values derived from the same data. If the sub- 
cluster had been continuously losing DM particles, we would 
expect an anomalously low M/L value for the subcluster — 
and in particular, a value lower than the main cluster's. From 
this agreement we can infer that the subcluster could not have 
lost more than / w 0.2-0.3 of its initial mass within r^- We 
should make two supporting comments here — first, the sub- 
cluster crossing time for its member galaxies is comparable 
to the time it took for the subcluster itself to cross the main 
cluster Thus, if the subcluster's gravitational well has been 
slowly diminishing due to DM collisions, its galaxies did not 
have time to evaporate, so this process can be ignored. Sec- 
ond, although the main cluster's DM particles are similarly 
affected by the subcluster impact, the effect on the main clus- 
ter's M/L value is smaller by at least their mass ratio (about a 
factor of 6). 

To escape beyond r = rn, the average subcluster particle 
needs a velocity Vesc — 1200 kms"'. The subcluster experi- 
ences a net loss of a particle in a collision if both vi > Vesc and 
V2 > v'esc, which occurs at scattering angles 6 given by (from 
eqs.|4ji; 

„2 \ 1/2 



< sin - < I 1 - 
vo 2 



esc 
„2 



(16) 



For our Vesc and vq, such angles exist. The probability of this 
loss, per collision, is 

J „ 2-71 f.m9d0 ,,2 



X = 



i,,^,2n sinOdO 

'(io) _ -j^ _ 2 esc 

2tt sin Ode ~ vl ' 



(17) 



where the upper integral is calculated over the interval of 6 
given by eq. il6i . During the subcluster's transit through the 
main cluster, each of its original particles has a probability to 
collide 

T-n,= -E„. (18) 

m 

Combining and ( I18> . the fraction of particles lost is 

21 



XT, 



m 



1-2 



v'o 



(19) 



where v'^^^. is a slightly higher escape velocity that allows for 
the subcluster's mass decline by up to a factor of 1 + / along 
the way; v'^^^ w Vesc(l + a/I +/)/2. We also ignore the subclus- 
ter's higher velocity at core passage, assuming the constant 
present velocity. (Both effects are small and our simplifica- 
tions are conservative.) Requiring that < f = 0.3, from 
eq. il9\ we obtain 



< 1 cm^ g ' . 



(20) 



The expression in the brackets in eq. jl9> is close to 1 for 
any small value of v'esc, thus our resulting constraint is only 
weakly dependent on the exact mass profile of the subcluster 
We also note that the subcluster's M/L ratios derived by C04 
represent lower bounds on the true values (because of possible 
projection and the particular technique employed to derive the 
subcluster mass), which makes our constraint conservative. 

For cross-sections such as (I20> and the subcluster's ob- 
served mass and radius, the mean scattering depth of the sub- 
cluster is Ts < 0.2. Thus our single-scattering assumption is 
reasonable as a first approximation. If one considers a pos- 
sibility that the escaping particle may expel another subclus- 
ter particle , the probability ( I19> is a conservative underesti- 
mate. In ^2.21 even higher Ts are allowed by our resulting 
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limit il5\ : for that method, single scattering is also a con- 
servative assumption. Inclusion of multiple scattering would 
lead to tighter limits, but requires detailed knowledge of the 
matter distribution inside the subcluster and is not warranted 
by the present data. 

3. DISCUSSION 

3.1. Improving the constraints 

From the above order-of-magnitude estimates, the most 
promising way of improving the limit on (or measuring) the 
DM collision cross-section using IE 0657-56 would be to re- 
fine the weak lensing mass map and limit (or detect) the mass 
loss from the subcluster core. Indeed, if it had lost a sig- 
nificant fraction of its DM particles, they should form a tail 
detectable in the dark matter map, not unlike the X-ray gas 
tail. This tail would, of course, be in addition to the matter 
tidally stripped from the subcluster and a gravitational wake 
discussed by Furlanetto & Loeb (2002). The main cluster's 
DM particles also would be scattered, probably resulting in 
additional detectable effects in the mass map. These consider- 
ations, along with the conservative assumptions that we made 
throughout to simplify calculations, suggest that the best way 
to proceed with the data interpretation is a detailed hydrody- 
namic simulation of this merger, such as that performed by 
Tormen, Moscardini, & Yoshida (2003) but with the inclu- 
sion of collisionless galaxies and coUisional DM with various 
cross-sections, followed by comparison of the results with the 
X-ray, optical and lensing data. 

3.2. Astrophysical context 

The nature of dark matter is still a great unsolved astro- 
physical problem. Laboratory searches for the DM candidates 
succeeded in showing that any interaction between the dark 
and baryonic matter is vanishingly small, with cross-sections 
many orders of magnitude lower than the values we consider 
in this work (e.g., Bernabei et al. 2003 and references therein). 
Despite their size, galaxy clusters have an average projected 
mass density of order 0.1-1 gcm"^, so they are no match for 
the laboratory experiments for constraining the DM-baryon 
interactions. However, clusters and galaxies may provide the 
best available laboratory for studying the DM self-interaction. 

While the common wisdom holds that DM is collision- 
less, a hypothesis of self-interacting dark matter (SIDM) 
with cross-sections of order 1-100 cm^g"' was most re- 
cently proposed by Spergel & Steinhardt (2000) to alleviate 
several apparent problems of the collisionless CDM model, 
such as the non-observation of predicted cuspy mass profiles 
in galaxies (e.g., Moore 1994; Flores & Primack 1994; cf. 
Navarro, Frenk, & White 1997; Moore et al. 1999b) and over- 
prediction of small sub-halos within the larger systems (e.g., 
Klypin et al. 1999; Moore et al. 1999a). Simulations and the- 
oretical studies (e.g., Dave et al. 2001; Ahn & Shapiro 2002) 
narrowed the range required to explain the galaxy profiles to 
a jm^ 0.5-5 cm g"' and pointed out that fluid-like DM with 
a/m^ 10"* cm^g"' is another possibility. 

Several observational constraints have been reported. 
Gnedin & Ostriker (2001; see also Hennawi & Ostriker 2002) 
pointed out that unless a/m < 0.3- 1 cm^g~', galactic halos 
inside clusters should evaporate on the Hubble timescale, be- 
cause of collisional heat conduction from the hot cluster DM 
particles into the cool halos. We note that our constraint in 
^2.3l uses a very similar idea, except that we have a cluster- 
sized halo subjected to a flow of particles from one direction 
instead of random bombardment. 



Furlanetto & Loeb (2002) proposed to use the different 
shapes of a gravitational wake that a subcluster halo moving 
through a bigger halo would create in the collisionless and 
fluid-like DM models. They argue that the X-ray image of the 
bright galaxy in the Fornax cluster akeady disfavors fluid-like 
DM. Their method is somewhat similar to our IE 0657-56 es- 
timates in that it also uses the motion of a halo within a bigger 
halo, but appears more observationally challenging. Follow- 
ing another suggestion of Furlanetto & Loeb (2002), Natara- 
jan et al. (2002) used observed sizes of the galactic halos in 
the A22 1 8 cluster (which should be truncated at different radii 
for different DM cross-sections) to obtain a/m< 40 cm-g~'. 
Upper limits in the 0.02- 10 cm^g"' range were derived by 
Hennawi & Ostriker (2002) from the absence of supermas- 
sive black holes in the centers of galaxies. 

Strong constraints are reported from the galaxy cluster 
cores. Yoshida et al. (2000) simulated evolution of a clus- 
ter for a/m = 10, 1, and 0.1 cm^g"' and obtained systemati- 
cally flatter radial mass profiles for higher a/m. Using Chan- 
dra X-ray data and the assumption of hydrostatic equilibrium, 
Arabadjis et al. (2002) derived a mass profile for the cluster 
MS 1358H-62 that is strongly centrally peaked (in fact, more 
peaked than even the collisionless simulations predict). From 
the comparison with Yoshida et al. (2000), they conclude that 
a/m < 0.1 cm-g"'. There are two potential difficulties with 
this stringent limit, however. Chandra revealed widespread 
gas sloshing in the cores of "relaxed" clusters (Markevitch, 
VikhHnin, & Forman 2002), including in MS 1358H-62, which 
does not lend support to the hydrostatic equilibrium assump- 
tion and the resulting X-ray mass estimates at relevant radii. 
More importantly, Yoshida et al. present time evolution of 
their simulated cluster profile for a/m = 10 cm^g"', which 
at different epochs covers a range of shapes comparable to the 
whole difference between their 0.1 and 10 cm^g"' simula- 
tions (see also Hennawi & Ostriker 2002). This suggests that 
for a high a/m, one expects a variety of cluster mass profiles 
— and indeed, while many real-world clusters exhibit peaked 
profiles, a large fraction has flat cores (e.g.. Coma). There- 
fore, a sample of cluster profiles is required for such stud- 
ies. Miralda-Escude (2002) used lensing data for the cluster 
MS 2137-23 to show that its central mass distribution is ellip- 
tical, and concluded that a/m < 0.02 cm^g"\ since otherwise 
DM collisions would erase ellipticity. Again, this interpreta- 
tion may not be unique, because ellipticity may be affected by 
line-of-sight projections. This ambiguity is illustrated by the 
fact that Sand et al. (2002), using lensing data for the same 
cluster but considering the flatness of the central mass distri- 
bution instead of its ellipticity, arrived at the opposite conclu- 
sion. 

While the above methods may eventually provide more sen- 
sitive constraints on a/ m than we have obtained, they require 
a statistical sample of clusters and cosmological simulations 
to reach solid conclusions. In comparison, the unique ge- 
ometry of the IE 0657-56 cluster merger allows us to see 
directly the (cumulative) results of single DM particle colli- 
sions, which makes the resulting limits quite robust. 

Finally, we note that our limit, a/m < 1 cm^g"', excludes 
most of the 0.5-5 cm-g~' interval proposed to explain the 
flat mass profiles in galaxies. Within the SIDM paradigm, 
the galaxy profiles and the tight cross-section limits coming 
from clusters can still be reconciled if the cross-section were 
velocity-dependent, so that it would be smaller on average 
in clusters than in galaxies (e.g., Firmani et al. 2000, 2001; 
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Hennawi & Ostriker 2002; Colin et al. 2002). However, it tions of this merging system, 
is difficult to justify this additional degree of freedom in the 
model until a nonzero cross-section is detected at any velocity. 



4. SUMMARY 

We have combined new X-ray, optical and weak lensing 
observations of the unique merging cluster IE 0657-56 to 
derive a simple, direct upper limit on the dark matter colli- 
sional cross-section, a jm < 1 cm^g"'. This is only an order- 
of-magnitude estimate; a more accurate, and quite possibly 
stronger, limit may be derived through hydrodynamic simula- 
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